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ABSTRACT 

Context. HIP 965 15 A is a double-lined spectroscopic binary included in the SACY catalog as a potential young star, and classified 
as an eclipsing binary in the ASAS catalog. It has a visual companion (CCDM 19371-5134 B, HIP 965 15 B) at 8'.'6. If bound to the 
primary, the optical and infrared colors of this wide companion are consistent with those of a white dwarf. 
Aims. The aim of this work is to characterize the system HIP 965 15 A&B by studying each of its components. 
Methods. We have analyzed spectroscopic and photometric observations of HIP 965 15 A and its visual companion, HIP 965 15 B. To 
confirm the system as a common proper motion pair, we have analyzed the astrometry of the components using high-angular resolu- 
tion infrared observations obtained within a time span of two years and archival astrometry. 

Results. The high-resolution optical spectrum of HIP965I5A has been used to derive a mass ratio, Mj/Mi, close to 0.9, 
with the components showing spectral types of Ml and M2. The optical light-curve of HIP 965 15 A shows periodic variations 
with Porb=2. 3456 days, confirming that HIP 965 15 A is an eclipsing binary, with preliminary parameters of !=89^0±ff2, and 
Mi=0.59±0.03 Mg and M2=0.54±0.03 Mg, for the primary and secondary, respectively, at an estimated distance of 42±3pc. This 
is a new eclipsing binary with component masses below 0.6 Mg. 

Multi-epoch observations of HIP 965 15 A&B show that the system is a common proper motion pair The optical spectrum of 
HIP 965 15 B is consistent with a pure helium atmosphere (DB) white dwarf. The comparison with evolutionary cooling sequence 
models provides reff,WD = 19, 126±I95 K, loggwD=8.08, MwD/Mg=0.6, and a distance of ~46pc. The estimated WD cooling age is 
~100Myr and the total age of the object (including the main-sequence phase) is ~400Myr. Finally, if HIP 96515 A&B are coeval, 
and assuming a common age of ~ 400 Myr, the comparison of the masses of the eclipsing binary members with evolutionary tracks 
shows that they are underestimated by ~I5% and 10%, for the primary and secondary, respectively. 

Key words. Stars: white dwarfs - Stars: binaries: eclipsing - Stars: individual (HIP 965 15) 

1. Introduction (d=43.9+9 pc, lvanLeeuwenll2007h Ml -ty pe star class i fied a s 

a double-line spectroscopic binary (SB2) bv lTorres et al.l ( l2006l) . 
In recent years, several infrared surveys have been focused The object has a visual companion at an angular separation of p~ 
on the detection of substellar companions around young 8^6 (CR 23, CCDM 19371-5134B, ' Worlev & Dougkiil fT^ 
and nearby stars usin^ high-contrast imaging techniques Dommaniit&Nyi I2QQ3 HIP 965 15 B hereafter). This visual 
_£e,g4Chauvin et al. 2003; L owrance et al. || 2005 |; [ Masciadri et alj binary has been detected in diff^erent infrared and optical sur- 
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iLafreniere et al. 2007). Young nearby stars are privileged ^^^^^ although its properties have never been studied in detail 



targets for these kind of studies due to their proximity (below 

100 pc) which allows to inspect very small separations, and their ^ ^''^^ ^tep to unveil the nature of the wide visual compan- 

youth, which enhances the contrast between the central star and i"" "IP 965 15 A, we compared our g -band photometry with 

the potential companion evolutionary tracks (' Baraff'e et al.l 120031; IChabrier et all I2000I) . 

In 2006 we started an infrared survey to detect s ubstellar ob- ^he comparison suggested that HIP 965 15 B may be of substel- 

jects around young stars from the SACY catalog ( iTorres et alJ "^^^^"^ *^ °f P^^^^y b^^^^^" 0.07- 1 Gyr and the 

Si l2006l) . Our sample included HIP 965 15, a nearby "^J^'^^.^ co-moving. However, the optical and infrared 

colors retrieved from archived observations are not consistent 

T , . , TZ 77 u with a substellar object but with a white dwarf companion. 

Send ojfprmt requests to: N. Huelamo ^ 

* Based on observations collected at the Paranal Observatory under Motivated by this preliminary result, we decided to perform a 

programs 77.C-0483(A) and 8 1 .C-0826( A). detailed study of HIP 965 1 5 . In this work, we have analyzed new 

** Marie Curie Outgoing International Fellow and archived observations of the triple system and we have de- 
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Table 1. Main properties of HIP96515 A. 



a(moo) s(moo) 



Spectral Parallax^ f^f^ f^^ 

Type [mas] [masyr'] [masyi->] 



19''37"08?7 -5r34'00:9 MP 22.79 91.07 -21.94 
±4.87 +4.93 +4.62 

Notes: ' Ivan Le euwen (2007): - iTorres et al. (2006) 

Table 2. Photometric data of the visual binary HIP 965 15 A&B 



Target 


V 

[mag J 


R 

[mag] 


J 

[mag] 


H 

[mag] 


[mag] 




HIP 965 15 A 


11.78' 


10.90' 


8.82- 


8.152 


7.962 






±0.04 


±0.03 


±0.03 


±0.04 


±0.03 




HIP 965 15 B 


13.0^ 


14.3" 


14.54- 


14.442 


14.522, 


14.94^ 






±0.3 


±0.18 


±0.51 


±0.31 


±0.08 



Notes: ' Torr es et al.l l|2006): 2 2MASS catalog ( Cutri et al. 2003); 
^ Visual Double stars in Hipparcos catalog dPonunanget & Nvsii2000D : 
" UCAC2 catalog; NACOA^LT data (this work). 



rived the main properties of each source. The paper is organized 
as follows: Section 2 describes the main observational proper- 
ties of HIP 965 15 A&B. In section 3, we analyze spectroscopic 
and photometric data of the primary star HIP 965 15 A, while the 
true nature of the visual component, HIP 965 15 B is unveiled in 
section 4. A discussion of the system is provided in Section 5, 
and the main results are summarized in Section 6. 

2. Observational properties of I-IIP96515 A&B 

The main stellar properties of HIP 965 15 A are displayed in 
Table [1] while its optical and infrared photometry are shown in 
Table |2] 

The analysis of the optical spectrum of HIP 965 15 A ob- 
tained during the SACY survey has allowed iTorres et al" l (l2006l) 
to classify it as an SB2, with almost equal-mass components (see 
TablelSj. From the cross-correlation function (CFF), we have es- 
timated a magnitude difference of A V- 0.55 mag between the 
binary components. 

HIP 965 15 is included in the ROSAT All-Sky Bright Star 
catalog ( Voges et al.ll999l) . The source was detected by ROSAT 
with a count rate of 0.172+0.036 counts s"'. We have de- 
rived the X-ray flux using the conversion factor definition by 
ISchmittetan(ll995h . The X-ray luminosity, assuming a distance 
of 44 pc (see Table 1), is Lx^(3 .61±0.75) lO^Vgs'' Using 
the bolometric correction (EC) by iKenvon & HartmannI (Il995h 
for an Ml -type star, and the V-mag listed in Table |2] we de- 
rive a log(Lx/i'boi) ratio of -3.07. If we assume that the two 
components of the spectroscopic binary (HIP 9651 5 Aa & Ab) 
are emitting in X-rays, so that L^Aa = LkAi, = L Jl, and us- 
ing the V-mag of each of the binary members (see Table 3) 
and the BC for a Ml and M2 star, respectively, we derive 
log(Lx/Lboi) ratios of -3.17 and -3.04 for Aa and Ab. These 
values are comparable to those fou nd in very active pre- main 
sequence (PMS) late-type stars (e.g. lBarrado v Navascues et alj 
[1999; Stelzer & Neuhauser 200 iV and in main -sequence (MS) 
late-ty pe stars in very close binaries (see e.g. iLopez-MoralesI 
l2007h . 

The age of HIP 965 15 A is not known. The lithium con- 
tent is normally used as a youth indicator in the case of late- 
type stars. However, it is not always pres ent in youi ig M- 
type stars since it is depleted too fast (e.g. iBasri et all Il996t 
lNeuhauseii[T997l: ITorres et"ani2008ah . In fact, none of the binary 
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Fig. 1. Hipparcos and ASAS V-band light-curves of 
HIP 965 15 A 



members of HIP 965 15 A show the lithium absorption doublet 
(/l=6707.8 A) in the optical spectrum. ITorres et alJ (l2008al) have 
studied the lithium content as a function of the V-I color in differ- 
ent young associations with ages between 5-70 Myr. The com- 
parison o f HIP 96515 A with similar M-type stars (with y-/~1.8, 
[Torres e t al. (2006)) in those associations allows to compute a 
lower limit to the age of > 70 Myr. By following t he temporal 
evolu tion of lithium in stellar models (eg. those bv iLandin et al.l 
l2006h for stars with the masses of HIP 965 15 A (Sect.©, Li 
should be exhausted by an age between 30 Myr (for the MLT 
parameter a-2) and 50 Myr (a= 1 ). 

Finally, HIP 9 65 15 B, is included in the AC catalog 
(IUrbanet al.lll998l) as a visual companion to HIP 965 15 A. The 
object has been detected in different infrared and optical surveys 
and multi-wavelength photometry is available in different pub- 
lic catalogs. We have included all the available data in Table |2] 
Note that, in most cases, the photometric measurements show 
large uncertainties. 

3. Characterization of I-IIP96515A: a l\11+l\12 
eclipsing binary 

In order to characterize the SB2 in more detail, we retrieved 
the phot ometric data of HIP 9 65 15 A from both the Hipparcos 
archive dPerrvman et al.|[r997l 96 points) and from the "All Sky 
Automated Survey", ASASPl ("Poima nski & Macieiewskil 120051 
818 points). The Hipparcos observations, in its own photometric 
system (Hip), were corrected with the relation Vj = Hip - 0.08, 
suggested for the V - / of this object (V - 7=1.85; [Torres et"aLl 
2006). 

The observations of the source are displayed in Figure [1] 
showing significant va riability. The object has been classified as 
an eclipsing binary bv iPaczviiski et al] (l2006h, using the ASAS 
data a lone. We have applied the method by i Lafler & KinmanI 
(119651) to search for any periodicity in the data. The best period, 
2.345 days, have been used to phase-fold both the ASAS and 
Hipparcos data, confirming the SB2 as an eclipsing system. In 
the next subsections, we describe our full analysis of the photo- 
metric data and derive the main properties of the eclipsing com- 
ponents, Aa and Ab. 



3. 1 . Photometric data: Period search 



The method by ILafler & KinmanI (ll965h was used to search for 
periods in the Hipparcos V-band observations of HIP 965 15 A. 
We searched for periods from 0.5 day to 63 days, by using a rel- 
ative period step of 6P/P=4xlQ^^ . Since M stars may flare, and 
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Table 3. Spectroscopic data of the eclipsing binary HIP 965 15 A 



Star 


V 


Spt. Type 


Vrad 


HaEW 


v.sin; 




[mag] 




[km s-'] 


[A] 


[kms-'l 


9651 5 Aa 


12.29 


MlVe 


84.1 


0.9 


14.8+1.5 


9651 5 Ab 


12.84 


M2:Ve 


-81.3 


2.0 


15.8±1.6 



HIP 965 15 is a strong X-ray emitter, we excluded the 4 bright- 
est measurements (from Hipparcos measurements) in this period 
search, in order to avoid these clearly non-periodic events, if they 
are present. Around 60 significant periods were found, but the 
most significant was f orb=(2.344451±0.000002)d. In the anal- 
ysis of the photometry with a version of the WD model (see 
Sect. l3.2b . we adjusted a phase shift to the folded observations, 
and obtained the ephemeris (HJD; the digits in parentheses affect 
the last digits) 



Min I = 2,448,493.80515(10) + 2.3444507(15)£, 



(1) 



which we use in our analysis. According to Eq, ([T) the secondary 
minimum occurs very close to phase 0.5, indicating an orbit with 
a small (if any) eccentricity, and in the lack of further evidence, 
we assume in the preliminary analysis, a circular orbit. We note 
that circular orbits are consistent with short orbital (and rota- 
tional) periods, like the one found in HIP96515A (e.g. |Z^ 
[1989). The V Ught-curves (Hipparcos and ASAS), folded with 
these ephemeris, are shown in Fig. |2] The stellar eclipses are 
clearly seen in both datasets. 

3.2. Preliminary analysis of the HIP96515A iight-curve 

We have analyzed the Hi pparcos optical light-curve with a 
version of the WD model dWilson & Devinnevlll97U Wilson 
[1979, 1993) extensively improved as described in IVaz et al. 
(|2007|) and references therein. Our version uses stellar atmo- 
spheres and models the radiated flux of both components us- 
ing the PHOENIX atmosphere models dAUard & Hauschildtl 
[T99llAllard etanfT997l:lHauschildt et al.lll997bllah . The (linear- 
law) limb-d^kening coefficients for both components have been 
taken from IClaretl (l2000h . and interpolated using a bi-linear 
scheme for the current values of log g and T^s at each iteration. 
The reflection albedos for both components have been fixed to 
a value of 0.5, which is appropriate for stars with convective 
envelopes. The gravity-brightening exponent, fi, has been com- 
puted using the local value of T^s for each point on the stellar 
surface and taking in to ac count mutual illumination following 
lAlencar& W3 (119971) and lAlencar et all (l999). Our procedure 
combines the WD modified code with several UNIX scripts and 
auxiliary programs to guarantee the consistency of the solutions 
in all steps. 

We have a single measurement of the radial velocity for each 
component at the phase 0.7715, according to Eq, ([TJ, by chance, 
very close to the quadrature (0.75 in a circular orbit), when the 
line separation is maximum. Since we do not have any estima- 
tion of the center of mass velocity of the system, the mass ra- 
tio is undefined. Therefore, to investigate the solution space, we 
have decided to make a grid of solutions for different values of 
the mass ratio, qi-MolMi), varying from 1.000 to 0.850 in steps 
of 0.025. By making this grid we can determine the velocity of 
the center of mass of the system for each q, and adjust the rest 
of parameters to reproduce the individual heUocentric velocities 
measured at phase 0.7715. 

We have assigned an arbitrary initial value of 89° to the sys- 
tem inclination, and we have used the WD model to find initial 
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Fig. 2. The phase-folded light-curves of HIP 965 15 A from 
ASAS {+ symbols) and Hipparcos (x symbols), consisting of 
818 and 96 observations, respectively, with our theoretical solu- 
tion for ^=0.9 obtained in Sect. |3.2| The Hipparcos observations 
were shifted by 0.6 mag, in order to separate the light curves. 
Note that ASAS observations were not corrected from the long 
term variation of the maxima, shown in Fig.[T] 

Table 4. Grid of solutions for 7 values of the mass ratio, 
q{^M2/Mi), for HIP 965 15 A. Note that Vi and V2, below are 
the theoretical velocities for the primary and the secondary at 
the orbital phase <p-Q. 11 15. All the solutions have the reduced 
the closest to 1 possible, for both data sets. 



<7(=M,/M,) 


1.000 


0.975 


0.950 


0.925 


0.900 


0.875 


0.850 


V, (km.s-') 


82.70 


81.65 


80.60 


79.48 


78.35 


77.19 


75.99 


V2 (kms"') 


-82.70 


-83.75 


-84.82 


-85.92 


-87.05 


-88.21 


-89.40 


(km ) 


1.40 


2.45 


3.52 


4.62 


5.75 


6.91 


8.10 


a:, (km s-') 


83.46 


82.41 


81.32 


80.21 


79.07 


77.90 


76.69 


K2 (kms-') 


83.46 


84.51 


85.60 


86.71 


87.85 


89.02 


90.23 


Ml (Me) 


0.56 


0.57 


0.58 


0.59 


0.59 


0.60 


0.61 


Ri (Ro) 


0.64 


0.65 


0.65 


0.64 


0.64 


0.64 


0.64 


loggi (cgs) 


4.57 


4.57 


4.58 


4.59 


4.60 


4.60 


4.61 


Ml (Mo) 


0.56 


0.56 


0.55 


0.54 


0.53 


0.53 


0.52 


Ri {Ro) 


0.53 


0.53 


0.53 


0.54 


0.55 


0.52 


0.52 


l0g«2 (cgs) 


4.74 


4.74 


4.73 


4.71 


4.69 


4.72 


4.72 


aiRe) 


7.73 


7.73 


7.73 


7.73 


7.73 


7.73 


7.73 


Am 


0.875 


0.893 


0.872 


0.832 


0.794 


0.891 


0.888 


in 


89.27 


89.45 


89.28 


89.08 


89.00 


89.51 


89.62 




42 


42 


42 


42 


42 


42 


41 



values to the gravitational pseudo-potentials (which define the 
size and form of the components) in order to reproduce the op- 
tical light-curve. Finally, the initial values of the effective tem- 
perature sjiave been adopted from the tabulated values listed in 
Poppdl ( 119801) . and Teff pri was kept fixed at the adopted value 
of 3,714±150K. In the lack of further information, the eclipsing 
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0.2 h 



0.0 0.? 0.4 0.6 0.8 1.0 

phose 

Fig. 3. The residuals O-C between the V light curves of 
HIP 96515 A from Hipparcos (x symbols) and ASAS (+ sym- 
bols) shown in Fig.|2l and our theoretical solution for q-Q.9 ob- 
tained in Sect. 13. 21 The standard error (1 observation) is shown 
in the upper left part of each panel. 



components were considered to be synchronized with the orbital 
motion, an assumption that may not be true (see below). 

Despite the presence of a visual companion to HIP 965 15 A 
at 8'.'6, we do not expect any contamination in the Hipparcos 
light-curve of the primary, since HIP 965 15 B is outside the 
Hipparcos field of view. On the other hand, the resolution of 
the ASAS setup is 14'.'2/pixel and the Ught from HD96515B 
is definitely included in their measurements. Therefore, we will 
assume a small amount of third light contribution in the anal- 
ysis of ASAS light curve. From the appearance of the folded 
light curve (using our ephemeris) and, in the lack of a well de- 
termined radial velocity curve, we have assumed a circular orbit 
for our analysis. 

Taking into account the above mentioned initial values and 
assumptions, we have modeled the orbital inclination, the grav- 
itational pseudo-potentials, the secondary effective temperature, 
a phase shift of the primary minimum and the primary internal 
luminosity for each point of the grid. The most interesting results 
from our modeling are displayed in Table |4] 

After convergence of the solutions of Table |4] we adjusted 
the amount of third light of the ^=0.900 solution, due to the 
evidence that the visual companion, HIP 965 15 B, contaminates 
the ASAS measurements. The third light parameter, in units 
of the eclipsing system light at quadrature, determined was 
X3=0.005±0.038 and the fit was essentially indistinguishable 
from the solution with X3=0. 

In order to reproduce the observed radial velocities, we have 
adjusted the semi-major orbital axis which, together with the 
orbital inclination and period, fixed the theoretical radial ve- 
locity amplitudes, Ki and K2 (in km s '). The apparent V- 
magnitude retrieved from SIMBAD is ll.lOmag. However, 
the Hipparcos measurement converted to Johnso n's V gives 
1 1.7 3 mag and our measurement is 11.78 mag CTorres et al.l 
I2OO6 ). value adopted for the distance estimation . By using bolo - 
metric corrections for both components from iPoppej (Il980h . 
we can derive the absolute dimensions and, finally, estimate the 
distance to be (42+3)pc. The interstellar absorption in such a 
short distance should be negligible. The parallax for the sys- 
tem is (22.79+4.87) mas measured by Hipparcos, which gives 
(44+9) pc, in very good agreement with our determination. 

All solutions of Table|4] are essentially indistinguishable 
from the curve fitting point of view, as can be seen by the rms 
of the residuals (cr, 1 obs) in magnitudes. The probable intrin- 
sic variability of one or both components makes it necessary 
to follow the system photometrically closely in time, preferably 



in different passbands and simultaneously with the spectroscopy 
for radial velocity measurements. Only with more precise light- 
curves in different photometric bands and with a precise radial 
velocity curve it will be possible to eliminate this degeneracy of 
the solutions in the mass ratio. For the purpose of this paper, we 
will assume that the components are in the main-sequence and 
make use of the mass-luminosity relation to fix q. 

Figure |2] shows our solution for ^=0.9 and Fig. [3] shows the 
corresponding residuals O-C. The steady and slow decreasing 
of the maxima present in ASAS observations was not corrected 
for in this analysi s. It is com mon that M dwarfs present varia- 
tions (e.g. .Rockenfeller et al.i,2006) mainly due to spots on their 
surfaces. However, we cannot say, with the data we analyzed, if 
the variability is due to one or to both eclipsing components. We 
searched for periodicities in the O-C residual s from 0.08 day to 
50 day s, using the phase dispersion method bv lLafler & KinmanI 
d 19651) and relative period steps of 6P/P-3.d-3, and found a sig- 
nificant period of (0.24965+0.00001) day (about6hrs), as can be 
seen in Fig.H] It is interesting to note that the variation pattern 
is present in both the Hipparcos and the ASAS data sets. The 
cadence of these observations (one or two observations a day, 
being common large intervals without observations) is not ade- 
quate for detecting these systematic variations in the O-C, with 
the annoying consequence of introducing a multitude of aliases 
in the analysis. It is, however, clear that the period of 0.24965 
day, if real, is not the only one: Fig.[T]gives indication for longer 
period variations and we, as a matter of fact, estimate at least an- 
other possible period in the O-C of (4,100+200) days, but with 
little significance level, due to the relatively short interval of the 
observations (~6,630 days). 

The period of 0.24965 day detected may be due either to pul- 
sation or to spots on the surface of one o r both of the eclips- 
ing components. Rockenfeller et al.l (|2006) list other M dwarfs 
with 6 hrs period. However, if the variation is due to spots on 
HIP 965 15 A, it would indicate that the spotted component(s) 
spins significantly (more than 9 times) faster that the orbital mo- 
tion (see e.g. CV Boo Torres et al. 2008b). This is probably not 
the case in HIP 965 15 A: using an orbital period of 2.345 days, 
and the radius of the binary components from Table 5, we es- 
timate rotational velocities of ~14km.s"' and ~12km.s"' for 
Aa and Ab, respectively. These values are comparable with the 
ones derived from the optical spectrum (Table 3, and assuming 
an inclination of 89 degrees), indicating that the system is syn- 
chronized. Therefore, the origin of the 6h-period remains unclear 
with the current data-set. 

By combining the results of TablelH and taking into account 
the maximum variations of the variables along the grid, we can 
derive the preliminary absolute dimensions for HIP 965 15 A, 
which are listed in Table|5] The masses of the binary com- 
ponents, HIP 965 15 Aa and Ab, are below 0.6 Mq. Hence, 
HIP 965 15 A is the fourteenth confirmed EB with component 
masses below 0.7 Mq (see e.g. IShkolnik et al.ll2008l) . which are 
important to calibrate theoretical evolutionary models. 

Finally, we note that the eclipsing system definitely de- 
serves more observations, both multi-band photometry with a 
cadence that allows determination of precise light curves, and 
spectroscopy, in order to obtain reliable radial velocity curves 
and spectral types of both components. If possible, HIP 965 15 B, 
should also be observed simultaneously, in order to give consis- 
tency to the analysis. 
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Fig. 4. The residuals O-C between the V light curves 
of HIP 965 15 A from ASAS (+ symbols) and Hipparcos 
(grey circles) shown in Fig. [3] folded with the period of 
(0.24965+0.00001) day. Note that the pattern is present in both 
Hipparcos and ASAS data sets (two cycles are plotted). 

Table 5. PreUminary physical parameters for HIP 965 15 Aa&Ab 



HIP 96515 



NACO/Ks 



Aa + Ab 



5" 



B 

o 



Fig. 5. VLT/NACO &-band image of HIP 965 15 A and its visual 
companion, HIP96515B, obtained in June 2008. The projected 
separation is of 378 AU at a distance of 44 pc. Note that the 
eclipsing binary (Aa & Ab) is not spatially resolved by NACO. 



Primary Secondary 



Absolute dimension: 

Mass (Mq) 0.59+0.03 0.54+0.03 

Radius (Ro) 0.64+0.01 0.55 + 0.03 

logg(c.g.s.) 4.59+0.03 4.69+0.05 
Photometric data: 



log reff(K) 3.571+0.018 3.555+0.019 

logL/Lo -1.15+0.07 -1.35+0.08 

Mm 7.63+0.18 8.13+0.20 

My 9.10+0.18 9.90+0.20 

Lsec/ipri 0.48+0.24 

"orb 

(days) 2.3456 

i (deg) 89.0+0.2 

Distance (pc) 42+3 



4. Characterization of I-IIP96515 B 

In the next subsections, we will focus on the characterization of 
HIP 965 15 B. First, we have investigated if the object forms a 
common proper motion pair with HIP 965 15 A. To do this, we 
have analyzed diffraction-limited infrared observations of the 
pair obtained in the course of our program to detect substellar 
object around SACY targets (Huelamo et al. 2008, in prep.). As 
a second step, we have obtained an optical spectrum to derive 
the main properties of the object. 

4.1. NACO/VLT observations ofHIP96515 

We observed HIP 965 15 twice with NAOS-CONICA (NACO), 
the Adaptive Optics facility at the Very Large Telescope (VLT), 
the nights of the 25th May 2006 and 16th June 2008. The 
target was observed with the visible wavefront sensor and 
the Ks-band filter. We used the S27 objective (nominal plate 
scale of ~0.027"/pixel) which provides a total field-of-view of 
~27"x27". The total exposure time was ~12 minutes on source. 

The data were obtained using a random jitter between the 
exposures to compute the sky emission. The i mages have been 
reduced with the Eclipse reduction package dDevillardlll997h 



following the standard steps: dark-subtraction, flat-field division 
and image shifting and stacking. The final image of HIP 965 15 
is displayed in Figure |5] Apart from HIP 965 15 A, we only de- 
tected the akeady known visual companion, HIP 965 15 B, at a 
projected separation of 8'.'6 (~378AU at 44 pc) from the pri- 
mary. Note that the eclipsing binary is not spatially resolved in 
the NACO image. 

The difference in magnitude between HIP96515A&B is 
6.98+0.07 mag. Taking into account the brightness of the pri- 
mary {Kfi-7.9mag, see Table 2), we can estimate a /Tj-mag of 
14.94+0.08 for the secondary. 

The NACO images have been used to derive the separation 
and position angle of the two objects in two different epochs. 
To derive accurate astrometry of HIP 965 15 B in the two im- 
ages, we have calibrated the plate scale and orientation of the 
infrar ed detector, CONICA, with the astrometric calibrator 0' 
Ori C dMcCaughrean & Stauffedl994l) . The results are displayed 
in Table |6l 

4.2. HIP96515 A & B: A common proper motion pair 

The astrometric accuracy of the NACO/VLT data is high enough 
to investigate if HIP 965 15 A&B are co-moving. For complete- 
ness, we have compared the NACO/VLT data with older (and 
less accurate) astromet ric meas urements from public surveys 
(e.p. lDommanget & Nv s 2000; C utri et al.l 2003). The visual bi- 
nary is include d in the Washington Double Star Catalog (WDSC, 
IWorle y & Douglass 1997) with three observations between 1910 
and 2000. Table |6] includes the angular separation and position 
angle of the two components at five different epochs. As seen, 
both remain almost constant within a time span of 98 years. 

Figure |6] (top panel) shows the difference in right ascension 
(RA) and declination (DEC) of the binary members measured 
in five epochs between 1910 and 2008. We have also overplot- 
ted the expected difference in RA and DEC of a background 
object taking into account the proper motion of the primary. If 
we focus on the NACO/VLT data only (Fig. |6] bottom panel), 
which provides the best astrometric accuracy, we can confirm 
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Table 6. Relative astrometry of HIP 965 15 A & B 




HIP96515 B 



2000 

1974 

2006 
1910 ' 2008 



2008 (IF BKG) 
2006 (IF BKG) 



1974 (IF BKG) 



1910 (IF BKG) 



-5000 -10000 

Aa (mas) 



16/6/2008 (IF BKG) - 



16/6/2008 
23/5/2006 



-6250 ^ 
-5950 



-6050 -6100 
Aa (mos) 



Fig. 6. Differences in right ascention and declination between 
HIP 965 15 A & B at different epochs (black circles). The ex- 
pected position of a background object (grey circles and labels) 
at each epoch is also displayed. Top: All the astrometric data 
from 1910 to 2008. Bottom: Only the NACO data, which pro- 
vides the best astrometric accuracy. The figure shows that the 
two objects are co-moving, that is, they are a common proper 
motion pair. 



that HIP 965 15 A&B are comoving , that is, they form a com- 
mon proper motion pair 



4.3. The Optical Spectrum of HIP96515B: observations and 
modeling 

To unveil the nature of HIP 965 15 B, we have obtained an opti- 
cal spectrum of the target with the ESO Multi-Mode Instrument 
(EMMI) on the New Technology Telescope (NTT) on October 
2007. The wavelength range goes from 380 nm to 700 nm, with 
a spectral resolution of R~l 100. The spectrum has been reduced 
with standard routines within IRAF, including bias subtraction 
and flat-fielding, and has been flux-calibrated with a spectropho- 
tometric standard observed during the night. 



Obs. Date Separation PA 
[arcsec] [deg] 



Reference 



2008-06-16 8.586±0.008 224.3±0.2 NACOA'LT, this paper 

2006-05-23 8.586±0.010 224.2±0.2 NACOA'LT, this paper 

2000-10-10 8.3 ±0.2 224.6 2MASS; CutrietaL (2003) 

1974.620 8.4 223 Dommanget c& Nvs (2000) 

1910.705 8.9±0.3 222 Worley & Douglass (1997) 



HIP 96.51. "5 B 



10 

X 


026.18- 
120.82 

387.92 
471 .47 
1517] - 

713.1.-;. 

921.9? 
015.57- 
047.73 

875.60 


> 8 
E 

b 

" -1 
2 


Iff \MfA ^ - 

^1 I'll 




4000 4500 5000 5500 6000 6500 


1 ,2 


I T,„ (K) = 19126 \ 


1 .0 


ll yklto log g (c.g.s) = 8.08 1 


« 0.8 


IllfV vJl Moss(Me) = 0.65 


E 

o 0.6 


1 ^ 11 1 ^^'''A Jikjaiu j 


0.4 






4000 4500 5000 5500 6000 6500 
Wavelength (A) 



Fig. 7. Top: NTT/EMMI Optical spectrum of HIP 965 15 B. We 
have marked all the identified He I lines. Bottom: Fit of the op- 
tical spectrum using pure helium model atmospheres. The best 
fit parameters are displayed at the top right corner 

Table 7. Best fit parameters for HIP 965 15 B 



Teff (K) 


19126±195 


log g (cgs) 


8.08+0.07 


Mass (Mq) 


0.648 


Radius (Rq) 


0.01212 


My (mag) 


10.89 


Mj (mag) 


11.33 


distance (pc) 


46 


Progenitor Mass (Mq) 


3.1+0.8 


WD Cooling age (Myr) 


106 ± 13 


MS lifetime (Myr) 


- 300 ' 


Total age (Myr) 


-400 



see text 



The optical spectrum of HIP 965 15 B is displayed in 
Figure [T] The continuum increases towards bluer wavelengths 
and the main detected spectral features are the He I absorption 
lines, typical of helium (DB) white dwarfs. 

We have compared the optical spectrum of HIP 965 15 B 
with atmospheric models of DB white dwarfs, in order to de- 
rive the physical parameters the object. Our grid of model at- 



Huelamo et al.: The triple system HIP 965 15 



7 



mospheres and syntheti c spectra for DB stars is described in 
iBeauchamp et alj (Il996h . These include the improved Stark pro- 
files of neutral helium of [Beauchamp et al. ( 1997). The model 
atmospheres used here assume a pure helium composition and 
cover a range of T^s = 10, 000 (1000) 16, 000 (2000) 30, 000 K 
and logg = 7.0 (0.5) 9.0. Our fitting technique relies on 
the nonlinear least-squares method of Levenberg-Marquardt 
jPres s et al. 1986), which is based on a steepest descent method. 
The model spectra (convolved with a Gaussian instrumental pro- 
file) and the optical spectrum are first normalized to a continuum 
set to unity. The calculation of is then carried out in terms 
of these normalized line profiles only. We adopt a pure helium 
composition based on the lack of any Ha absorption feature; ad- 
ditional models with mixed helium and hydrogen compositions 
allow us to set a limit of A^(H)/A^(He) < 10"^ 

As usual, there is a cool and a hot solution on each side of 
the temperature at which the line strengths reach their maximum, 
but the hot solution near Teff =25,000 K can easily be ruled out 
based on the slope of the observed energy distribution. Our final 
solution, reff=19,130+190 K and logg=8.08±0.07, is displayed 
in Figure|2] These atmospheric parameters can be converted into 
a mass of M=O.65+O.O4M0 using evolutionary models with C/O 
cores, q(}ie) = log MHe/M* = 10"^, and ^(H)=10"'°, which are 
representative of helium-atmosphere white dwarffl 

By using the best fit parameters of the model we can de- 
rive the absolute magnitudes in different photometric bands 
(see Table |7]i. Using Mj and the 2MASS 7-band brightness 
of HIP 965 15 B, which is the measurement with the smallest 
photometric error, we estimate a the distance to the source 
of (i=45.8pc, which is in good agreement with both the 
Hipparcos parallax and the distance derived from the analysis 
of HIP 965 15 A. 

The evolutionary sequences described above, as well as em- 
pirical relations (see e.g. Fontaine et al. 2001), allow to esti- 
mate the mass of the white dwarf progenitor and its cooling age. 
These data are provided in Table |7] We obtain two different es- 
timates of the mass of the white dwarf progenitor. First, we use 
the em pirical initial-final mass relationship (IFMR) from lWoodI 
(|1992|) . which is based on the study of the wh ite dwarf luminos- 
ity function and m ass distribution (see also lLeggett et al.ll 19981 
iFontaine etalj2001i) . This yields a value of Mms = 3.8+0.5 Mp 
Altern atively, we use the IFMR determined by iFerrario et alj 
(12005') based on white dwarfs observed in open clusters (see 
in particular their Fig. 1). This second relation yields a signif- 
icantly smaller final mass of Mms - 2.5 + 0.5 Mq. However, as 
discussed by Ferrario et al., the IMFR of most clusters appears 
to have an intrinsic spread (i.e, there is a range of initial masses 
that gives rise to a given final mass). Hence a value of Mms as 
high as 3.8 Mq cannot be necessarily ruled out. In the following, 
we thus adopt a mean value of Mms = 3.1 + 0.8 M© to allow for 
the possible mass range of the progenitor. 

We have derived the main-sequence lifetime of HIP 965 15 B 
using two sets of evolutionary tracks for a 3.1+0.8Mq star 
First, we have used the ATON 2.4 models, which are a 
modified version of the models by D'Antona & Mazzitelli 
dP' Antona & Mazzitellil 11994 119971) . and are fully described 
in^andin et al. (2006). According to these evolutionary models, 
a 3.1 +0.8 Mq star spends 280;^j25Myr on the main-sequence, 
i.e., from the zero-age main-sequence (ZAMS) t o the termina l 
age main-sequence (TAMS). On the other hand, IClaretl (|2004 
models (Z=0.02) predict a span of ~317Myr for the MS time 
of a 3.1 Mq. Assuming a MS lifetime of ~300Myr and a cooling 



^ see |http://www.astro.umontreal.ca/~bergeron/CoolingModels| 



time of ~ 1 00 Myr (see Table|7]l, the total age of the system would 
be ~400 Myr, being the errors dominated by the uncertainty in 
the mass. 



5. Discussion 

Th e comparison of HIP 965 15 Aa & Ab with evolutionary tracks 
bv lBaraffe et al.l (Il998l) provides an age of ~50 Myr for masses 
of 0.5 Mq (with alphas 1), and places the binary on the PMS (see 
Fig. ^. The PMS nature of HIP 965 15 A is also suggested by 
tracks generated with the ATON 2.4 stellar evolutionary code, 
where we adopted the solar chemical composition (Z=0.0175 
and Y=0.27) and alpha, a, equal to 1. The tracks in Fig.[8]gener- 
ated with ATON2.4 code cover the PMS, MS and eai'lier stages 
of post main-sequence phases. The best age, according to ATON 
2.4, is ~60Myr (Fig. [8]), although with large uncertainties. 

On the other hand, our astrometric study shows that the 
white dwarf and the EB are co-moving. If we assume that 
HIP 965 15 A&B are coeval, the age of the triple system should 
be ~400Myr, which is significantly older than the age provided 
by the PMS tracks for HIP 965 15 A, although within the obser- 
vational uncertainties of the primary star (see Fig. [8]). 

To shed light on the evolutionary status of HIP 965 15 A, we 
have studied the kinematical properties of the object. Using an 
average distance of 45 pc (coming from the EB and W D inde- 
pendent studies presented here), the RV derived by Torres et all 
(2006), and the proper motions measured by Hipparcos (see 
Table 1), we can derive the UVW components of the Galactic 
space velocity vector: UVW=-6.Q, -0.5, -21.6kms"^ We have 
compared them with the UVW vectors of y oung moving groups 
in the sol ar neighborhood i ncluded in I Zuckerman & Son3 
(12004 and lTorresetaP (l2008ah . We conclude that HIP 965 15 A 
does not share the kinematical properties of any of these associ- 
ations. 

Assuming an age of 400 Myr for the triple system, a possible 
explanation to the age of 60 Myr derived with PMS evolution- 
ary tracks is the different evolution of EBs and very magnet- 
ically active stars in comparison with non-active objects. EBs 
are normally fast rotators (Prot < 3 days) and they exhibit strong 
X-ray emission related with the presence of magnetic fields 
dL opez-Morale si l2007h . As explained by e.g. IChabrier et al.l 
(20017 rapid rotation and/or magnetic fields reduce the effi- 
ciency of large-scale thermal convection in their interior leading 
to less efficient heat transport. The reduction of the stars radiat- 
ing surface due to the presence of cold spots can yield to smaller 
Teff and larger radii, thus dep arting significantly from the models 
predictions (see e.g. lTorre"s & Ribas 2002). 

To investigate this in more detail, we have studied the Radius 
vs. Mass relation for HIP 965 15 A and compared it with th e pre- 
dictions from evolutionary models by iBarafFe et al.l ( Il998l) . The 
adopted models correspond to a solar metallicity of [Fe/H]=0, 
and a mixing length of q'=1. The results are displayed in Fig.|9] 
which includes all the low-mass eclip sing binaries confirm ed so 
far The data has been adopted from IShkolnik et al.l (12008.) and 
references therein. 

If we assume an age of 400 Myr for the triple system, the 
radii of HIP96515 Aa&Ab are -15% and 10% larger, respec- 
tively, than predicted by the evolutionary tracks (see Fig. |9]l. 
These discrepancies have been previously observed in similar 
eclipsing binaries and have been attributed to the e ffect of rapid 
rotation, magnetic activity and/or metalhcity (e. g. iBerger et al] 
120061: IChabrier et alJl2007HL6pez-Moralesl2007h . 
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log Teff 



Fig. 8. Stellar Luminosity (top panel) and surface gravity (bot- 
tom panel) vs. effective temperature. HIP 965 15 Aa and Ab 
are represented by filled circles. The grey and blac k solid 
line s represent evolutionary tracks by iBaraffe et alJ (Il998h 
and lLandin et alj (f2006), respectively, for 0.6 and 0.5 M© stellar 
masses. The dotted lines represent the uncertainties in the stellar 
masses for the ATON models. The comparison of HIP 965 15 A 
with the two sets of evolutionary tracks provides an age of 
~60 Myr, that is, places the eclipsing binary members on the pre- 
main sequence, although with large uncertainties. 

6. Results and conclusions 

In this paper we have studied in detail HIP 965 15, a triple system 
composed by a double-lined spectroscopic binary (SB2) and a 
co-moving white dwarf at a projected separation of 8'.'6. We have 
analyzed both new and archived observations of the system. Our 
main results can be summarized as follows: 

- The analysis of multi-epoch optical photometry of 
HIP 965 15 A has revealed that this SB2 is an eclipsing 
binary with almost equal-mass components: M^a =0.59 M© 
and Ma/,=0.54Mo. In fact, this is the fourteenth confirmed 
eclipsing binary with component masses below 0.7 M©, 
which are of extreme importance to calibrate theoretical 
evolutionary tracks. 

- Archival astrometry and new NACO/VLT data show that 
HIP 965 15 A&B is a common proper motion pair. 

- The optical spectrum of HIP 965 15 B is consistent with a 
hot DB white dwarf. We have modeled it using evolution- 
ary cooling sequence models, deriving a temperature, sur- 
face gravity and mass of Teff (K)=19126, \ogg = 8.08, 



T 



0.8- 




0.0 0.2 0.4 0.6 0.8 

Mass (Mb) 



Fig. 9. Mass-Luminosity relation for all known eclipsing bina- 
ries with masses M < 0.7 Mq (see Shkolnik et al. 2008, and ref- 
erences therein). HIP 965 15 Aa and Ab are represented by filled 
circles. T he dotted, dashed and solid lines represent evolutionary 
tracks bv lBaraffe et"an (Il998l) for 50, 100 and 500 Myr, respec- 
tively. 

and M=0.65Mo. Empirical relations provide a progenitor 
mass of 3.I+O.8M0, and a white dwarf cooling age of 
100+13 Myr. The MS lifetime of a 3.1 Mq star, as derived 
from evolutionary models, is ~300Myr, resulting in a total 
age of ~40 Myr for the white dwarf. 

- A dire ct comparison of the EB members with evolutionary 
tracks (iBaraffe et al.lll998l: llTandin et al.ll2006h places them 
on the PMS with an age of ~60Myr, although with large 
uncertainties. We note that HIP96515 A does not share the 
the kinematical properties of known young moving groups 
in the solar vicinity. 

- If HIP 965 15 A&B are coeval, the system must be older 
than 300 Myr. The discrepancy between this age and the 
60 Myr provided by PMS tracks might be explained by the 
different evolution of EBs, which normally show fast rota- 
tion and strong magnetic fields, compared to non-magnetic 
stars. Since HIP 965 15 A shows fast rotation and strong X- 
ray emission, it can depart from the predictions of theoretical 
models, which place the binary on the PMS. 

- If we assume an age of ~400 Myr for the triple system, the 
comparison of HIP 965 15 Aa&Ab with evolutionary tracks 
by Baraffe et al. ( 1998) shows that the models underestimate 
the stellar radii by 15% and 10%, for the primary and the 
secondary, respectively. This behavior has been previously 
observed in similar low-mass EBs. 

Finally, we note that additional RV observations are needed 
to better constraint the orbit and physical parameters of the 
eclipsing binary. 
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